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ELECTRICAL NETWORKS AND STEPHENSON’S CONJECTURE 


SA’AR HERSONSKY 

We dedicate this article to the memory of Bill Thurston. 


Abstract. In this paper, we consider a planar annulus, i.e., a bounded, 2-connected, Jor¬ 
dan domain, endowed with a sequence of triangulations exhausting it. We then construct 
a corresponding sequence of maps which converge uniformly on compact subsets of the do¬ 
main, to a conformal homeomorphism onto the interior of a Euclidean annulus bounded 
by two concentric circles. As an application, we will affirm a conjecture raised by Ken 
Stephenson in the 90’s which predicts that the Riemann mapping can be approximated by 
a sequence of electrical networks. 


0. Introduction 

0.1. Riemann’s Mapping Theorem and Thurston’s disk packing scheme. The Rie¬ 
mann Mapping Theorem asserts that any simply connected planar domain which is not the 
whole plane, can be mapped homeomorphically by a conformal mapping onto the open unit 
disk. That is, the domains are conformally equivalent. After a suitable normalization, this 
mapping is called the Riemann mapping and it is desirable to have a concrete approximation 
of it. In [IH], Rodin and Sullivan proved Thurston’s celebrated conjecture [56] asserting that 
a scheme based on the Koebe-Andreev-Thurston disk packing theorem (cf. [H |21 SSI EZj) 
converges to the Riemann mapping. 

In order to formulate Thurston’s conjecture which inspired Stephenson’s conjecture, we 
need to recall a few dehnitions. Let P be a disk packing in the complex plane C. An interstice 
is a connected component of the complement of P; and one whose closure intersects only 
three disks in P is called a triangular interstice. We will let supp(P) denote the union of 
the disks in P and all its bounded interstices. The disks of P that intersect the boundary of 
its support are called boundary disks. Two hnite disk packings P and P in C will be called 
isomorphic if there exists an orientation preserving homeomorphism 0 : supp(P) —?> supp(P) 
such that 0(P) = P. It is evident that such an isomorphism h induces a bijection between 
the disks of P and the disks of P and that an isomorphism of packing is a combinatorial 
notion. 

Let C C be a bounded simply connected domain, and let po be an interior point in it. 
For each positive integer n, let P"^ be a disk packing in in which all bounded interstices 
are triangular. Assume that there is a sequence of positive numbers which converges to 
zero, such that: i) the radius of every disk in P"’ is smaller than 6n, and ii) every boundary 
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disk in P” is at distance at most 5n from dVL. Finally, let Pq be a selected disk in P” which 
is closest to po or contains it. 

The Disk Packing Theorem (Koebe-Andreev-Thurston) implies that there exists an iso¬ 
morphic packing P” in the closed unit disk D with all of its boundary disks tangent to 
the unit circle Furthermore, if the given graph is isomorphic to the 1-skeleton of a 
triangulation of the Riemann sphere, then the packing is unique up to applying a Mobius 
transformation. Let 

(0.1) /„ : supp(P") ^ supp(P") 

be an isomorphism of P” and P"^. Furthermore, normalize P” by a sequence of Mobius 
transformations preserving U so that Pq, the disk corresponding to Pq, is centered at the 
origin. Thurston conjectured that if the packings P"' are chosen to be sub packings of scaled 
copies of the inhnite hexagonal disk packing of C, then the sequence of piecewise affine maps 
(i.e., simplicial) /„ converges uniformly on compact subsets of D to the Riemann mapping 
from D to D. 

Rodin and Sullivan |18] proved Thurston’s Conjecture by hrst showing that the maps 
fn are iF-quasiconformal, for some fixed K. Hence, there exists a subsequence which will 
converge to a limit function / which must also be iF-quasiconformal. Rodin and Sullivan 
further showed that / must be 1-quasiconformal, and therefore, / is in fact conformal. He 
and Schramm [5^1 Theorem 1.1] developed profound techniques which avoid the machinery 
of quasiconformal mapping that is heavily used in Rodin-Sullivan’s proof. Up to date, their 
theorem and advances [3l] in the simply connected case, is the most advanced (see also their 
related work in [35]). 

More recently, Chow and Luo [TS] found applications of circle mappings to the study of 
discrete Ricci flow on surfaces; see also the work of Glickenstein m for related study. There 
are also applications of circle packings to algorithmic computer vision and computational 
conformal geometry due to Gu, Luo and Yau, Gu, Zeng, Luo and Yau, and Sass, Stephenson 
and Brock (cf. [3T1[32] and |19] as examples and further advances). 

0.2. Electrical networks and Stephenson’s conjecture. In his attempts to prove uni- 
formization, Riemann suggested considering a planar annulus as made of a uniform conduct¬ 
ing metal plate. When one applies voltage to the plate, keeping one boundary component 
at a hxed voltage k and the other at the voltage 0, electrical current will flow through the 
annulus. The equipotential sets form a family of disjoint simple closed curves foliating the 
annulus and separating the boundary curves. The current flow sets consist of simple disjoint 
arcs connecting the boundary components, and they as well foliate the annulus. Together, 
the two families provide curved “rectangular” coordinates on the annulus that can be used 
to turn it into a right circular cylinder, or into a (conformally equivalent) circular concentric 
annulus. 

An electrical circuit or network is a collection of nodes and connecting wires. For in¬ 
stance, a disk packing of a fixed planar domain induces such a network where each center 
of a disk corresponds to a node and a wire connects each pair of nodes corresponding to 
tangent boundaries. It is therefore reasonable to conjecture that if the domain is made of 
thin conducting material then its electrical behavior can be approximated by a sequence of 
networks that approximates its shape. 
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Stephenson’s conjecture from the 90’s is concerned with constructing such an approxima¬ 
tion: 

Conjecture 0.2 (Stephenson [SI]). Given a sequence of networks approximating a simply- 
connected, hounded, Jordan domain arising, for instance, from a sequence of disk packing, 
choose conductance constants along the edges (for each network) according to Equation / (i.gl) . 
Then the sequence of discrete potentials and currents will converge to the ones induced by 
the Riemann mapping. 

We have phrased this conjecture in the more recent formulation of fll.51) (see Section [1] for 
the details). In fact, a similar conjecture can as well be formulated for any domain that can 
be approximated (in a sense that we will dehne in Section 13.ip by a sequence of (rehned) 
quasi-uniform triangulations (see Dehnition IA.7p . 

In Theorem 13.121 we will formalize and affirm Stephenson’s conjecture in the case of a 
polygonal annulus by methods that are different from the ones used in his paper or those 
mentioned in Section 10.11 In particular, we will explore a large class of networks for which 
it holds. As one application, we will affirm Stephenson’s conjecture in its original form. 

0.3. The themes of this paper. There is a classical and elaborate theory of conformal 
uniformization for domains in the Riemann sphere that are bounded by non-degenerate 
Jordan curves. Let be such a domain which is also hnitely connected. Koebe proved [12] 
that U is conformally homeomorphic to some domain fl* whose boundary components are 
circles. Such a domain is called a circle domain. Furthermore, U* is unique up to Mobius 
transformat ions. 

Discrete uniformization schemes have traditionally been the hrst step in constructing a 
sequence of approximations to a conformal map from the given domain (more on this in 
Section EH). There is much current interest and effort by (for example) Cannon, Floyd and 
Parry to provide sufficient conditions under which, discrete schemes based on the discrete 
extremal length method, will converge to a conformal map in the cases of triangulated annulus 
or a quadrilateral (see for instance [Slilolllllill]). Of much current interest is the universality 
of the critical Ising and other lattice models where discrete complex analysis on graphs played 
crucial role (see for instance [I5l[23]). However, it was shown by Schramm [501 page 117] that 
if one attempts to use the combinatorics of (for instance) the hexagonal lattice alone, square 
tilings (as constructed by Schramm’s method) cannot be used as discrete approximations for 
the Riemann mapping. 

In this paper, stemming from our work in [SG] [S^ |38], we will prove that a certain discrete 
scheme yields convergence of the mappings described below to a canonical conformal mapping 
from a given polygonal, planar, annulus, onto the interior of a Euclidean annulus bounded 
by two concentric circles. 

Specihcally, the underlying idea of this paper is rooted in a foundational feature of two 
dimensional conformal maps. If / : D —)■ C is a conformal map, then the Cauchy-Riemann 
equations imply that 3f?(/) and A(/) are harmonic functions, and that A(/) is the harmonic 
conjugate of 3?(/). For instance, when {r,9) denote polar coordinates in the plane, we have 
that u{r,6) = logr and v{r,6) = 6 (when 6 is single valued) are harmonic functions, and 
v{r, 9) is the harmonic conjugate of u{r, 9). Indeed, in [38], starting with the solution g oi a. 
special discrete Dirichlet boundary value problem dehned on a polygonal annulus we dehned 
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two new functions on The first function, g*, which is actually dehned for the annulus 
minus a sh|3, is obtained by integrating the discrete normal derivative of g along its level 
curves (see [3H1 Dehnition 2.6]). The second function, h, which depends on g* is obtained 
via a solution of a discrete Dirichelt-Neumann boundary value problem on the domain of g* 
(see [3H1 Dehnition 3.1]). 

Hence, in the case of a planar triangulated annulus, and for any sequence of triangulations 
{Tn\ of it, we have combinatorial approximating candidates for u{r,9) (the gns), and two 
combinatorial candidates for v{r,9) (the g^^s or the h^s). However, in this paper, due to 
the current state of the art in the Loo convergence results which we will apply, we will work 
with the pair {g,g*) where g* (see Dehnition I2.32p is a signihcant modihcation of g*. In 
Theorem 13.121 we will show that, for a suitable sequence of quasi-uniform triangulations 
endowed with Stephenson’s conductance constants dehned according to Equation fll.2p . the 
affinely extended maps of {gn, will converge uniformly on compact subsets of a given 
annulus, to the real and imaginary parts of a conformal uniformizing map of the annulus, 
respectively. To this end, we will employ recent advances from the theory of the hnite volume 
method, techniques from discrete potential theory, and classical results form the theory of 
functions of one complex variable and partial diherential equations. 

In order to put the needed advances over previous applications of the hnite volume method 
in context, let us briehy recall an inspiring work by Dubejko [22]. Set w to be the solution 
of the Dirichlet problem Aw = f for x E Ll, and w = (f) for x G cID, where D is a simply- 
connected, bounded, Jordan domain with boundary, where / G L^(D) and (j) G C°(D). 
By applying techniques from the hnite volume method, Dubejko proved that w can be 
approximated (in various norms) by a sequence of solutions of discrete Dirichlet boundary 
value problems. These solutions belong to a certain Sobolev space and are constructed via 
a sequence of triangulations (with special properties) that gets hner while exhausting D 
from the inside. Dubejko’s work, which utilized Stephenson’s conductance constants, is not 
sufficient for constructing approximations of conformal maps from Jordan domains. In fact, 
already in the simply connected case his techniques are not sufficient. This is due to the 
following reasons: First, his methods can be applied only under the assumption that the 
boundary of D is second, Dubejko did not address the problem of approximating a 
harmonic conjugate; third, Dubejko utilized Riemann’s mapping theorem in his proof. 

0.4. Organization of the paper. In Section [H we describe the conductance constants 
suggested by Stephenson, and express these in the way applied in Theorem 13.121 the main 
theorem of this paper. 

In Section [2l we present three novel dehnitions. We dehne the class of discrete asymptotic 
harmonic functions. Intuitively, a function in this class is almost harmonic on a scale deter¬ 
mined by the mesh of the triangulation. This class will be employed in the approximation 
process described in the main theorem of this paper. The flux following path in a given 
triangulation is associated to the amount of discrete flux crossing a given path in the one 
skeleton of the Voronoi cells of the triangulation. Finally, utilizing a flux following path, we 
dehne a conjugate function, g*, of a discrete harmonic or a discrete asymptotic harmonic 
function. 

slit is a chosen simple path in the annulus connecting the two boundary components on which g is 
monotone increasing. 
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In Section 131 we provide the proof of Theorem 13.12[ where the convergence of our proposed 
discrete scheme is shown to converge in the case of a polygonal annulus. We are then con¬ 
cerned with the approximation of the uniformization of an annulus with continuous Jordan 
boundary. Corollary 13.381 demonstrates that Theorem 13.121 coupled with a generalization of 
a compactness theorem due to Koebe and a diagonalization process, allow the weakening of 
the boundary regularity assumption of Theorem 13.121 from polygonal to continuous. 

Section |4] is devoted to the proof of Theorem 14.21 where we address the approximation 
of the uniformization of a bounded, simply-connected Jordan domain, the setting in which 
Conjecture 10.21 was hrst stated. The idea is to present the punctured domain as an increasing 
sequence of annuli. Thus, one can apply Theorem 13.121 to each annulus in the sequence. 
The existence of a converging subsequence of the maps obtained in each step to a bounded, 
conformal, univalent map follows the same rationale as in Corollarv l3.381 and we can therefore 
restrict attention to the case that the boundary of the domain is polygonal. The final 
ingredient in the proof is the Riemann’s removable singularity theorem. 

With the aim of making this paper self-contained to a broad audience, it contains an 
Appendix. In Appendix lA.ll and in Appendix IA.21 we collect a few important notations, 
dehnitions and theorems from the finite volume method that are applied in this paper. The 
reader who is familiar with this method, can skip these sections. However, Theorem IA.231 
which is quoted from [T3] is essential for the Loo convergence analysis results of this paper. 
In Appendix IA.31 we explain how Stephenson’s conductance constants naturally appear in 
the setting of the finite volume element method. 

Acknowledgement. We are indebted to Gerard Besson, Danny Calegari, Gilles Gourtois, 
Benson Farb, Bill Floyd, John Hubbard, Rich Schwartz, Ted Shifrin and Robert Varley for 
their interest in this research, patient listening and valuable discussions. Ridgway Scott, 
Thierry Gallouet and A1 Schwatz are heartily thanked for graciously sharing with the author 
their insights regarding the subtle analysis involved in numerical methods for convex and 
non-convex planar domains. 

1. Electrical networks induced by disk packings and Stephenson’s 

CONDUCTANCES 

Let us recall a few definitions and some notation from Ell [231 El] and [33] in order to 
express the conductance constants suggested by Stephenson. Let P be a euclidean disk 
packing of a domain D for a complex K, i.e., the contact graph of P is isomorphic to K. 
For an interior edge {u,v) e K, consider the tangent circles Cy,Cu as depicted in the hgure 
below. Let Cx, Cy be their common neighboring circles. 

The radical center, Wx, of the triple {cy, c„, Cx} of circles will denote the center of the circle 
that is orthogonal to c„, Cy and Cx and let Wy be the radical center of the triple {cy, c„, Cy}. 
Let Zu,Zy be the centers of Cu,Cy, respectively. Finally, for a vertex v, let Ry denote the 
radius of the circle Zy. 

Stephenson’s eonduetance of an edge is defined by (see also Definition lA.lOl and Equa¬ 
tion (IA.28p L 

'^x '^y I 


( 1 . 2 ) 


c(e) = c{u, v) 
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Figure 1.1. Constructing an edge conductance in a circle packing. 


It is illuminating to give a probabilistic interpretation to this quantity. Stephenson’s main 
idea was to chase angle changes at the centers of the circles, as radii changes while maintaining 
(new) disk packing. Given a euclidean circle packing, the effect of a small increase in the 
radius of one circle, say R^, is that the angle sum 0^ decreases, while the angle sums 0^^. 
at the neighboring vertices {^ 1 ,^ 2 , • • • ,Vk} increase. Some of the angle “distributed” by v 
arrives at Vj and must be passed along in order to keep a packing at Vj. Hence, has to 
be adjusted and we need to keep track of the angle changes of its neighbors, and so forth. 

In Euclidean geometry, the angles of any triangle add up to tt, so angles in this process will 
never get lost. In other words, the total angle leaving one vertex must be divided into portions 
and then distributed as angles arriving to its neighbors. This movement can be expressed 
as a Markov process, where the transition probability from v to vj, is the proportion of a 
decrease in angle at v that becomes an increase in angle at Vj. In this Markov process, the 
random walkers are the quantities of angles moving from one vertex to another. Thus, for a 
specihc neighbor u = Vj, the amount of angle arriving at 'ijju is given by as a result, the 
transition probability from n to m as described above is given by 


(1.3) 


p{v,u) 


dipu 

dRv 

E k dipvj 

i=l dRZ 


Also, for a vertex v G K, we let 


(1.4) 


p{v,u) 


c{v, u) 


It is remarkable that in 2005 (see [55l Section 18.5]) Stephenson showed that equality 
holds between these two Markov transitions, that is. 


( 1 . 5 ) 


p{v, u) = p{v, u), M ~ n. 
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2. Smooth harmonic conjugate functions and their combinatorial 

COUNTERPARTS. 

This section entails several key definitions and constructions. In the first subsection, we 
collect a few classical PDE existence results that go back to Poincare and Lesbegue. In 
the second subsection, we will assume that ^ is a fixed, planar, polygonal annulus endowed 
with a triangulation T. After recalling the definitions of the combinatorial laplacian and 
the normal derivative, we will turn to define the class of discrete, asymptotically harmonic 
functions. This class will enable us to define a conjugate function to a function which is 
either harmonic or asymptotically harmonic of some order (see Definition I2.18p . 


2.1. Strong solutions of the Laplace equation and smooth harmonic conjugate 
functions. Let D be a planar domain and assume that u G (7^(12) fl C°(r2) is the strong 
solution of the Dirichlet boundary value problem for the Laplace equation 


( 2 . 1 ) 


Am = 0 in D, 
u = h on dfl, 


where h is the trace of h G It is well known that a strong solution is a also a weak 

solution (see Appendix |A.2| or m Section 3.2]) for the details and the definition of 

The study of the existence of strong solutions of Dirichlet boundary value type problems 
has an interesting history. Poincare introduced the notion of barriers, and their importance 
was further recognized later by Lebesuge. A function w G A‘^(D) is called super-harmonic 
in D, if for any subdomain D' C D, and any harmonic function m in D', 

(2.2) w > u, in Dh 

Let ^ be a point in dQ, then a (7^(12) function w = is called a barrier at ^ relative to 
12, if w is super-harmonic in 12, it approaches 0 at and outside of any sphere about it 
has a positive lower bound in D. Two profound consequences of the existence of a barrier 
are the following. 


Theorem 2.3 ([HI Theorem III, page 327]). A necessary and sujficient condition that the 
Dirichlet problem for If is solvable for arbitrary assigned continuous boundary value values, 
is that a barrier for 12 exists at every point in dQ. 


Theorem 2.4 f [Hl Theorem IV, page 328]). Let u be the solution of the Dirichlet problem 
for 12 with cf G C^{dLl) boundary values. If ^ E dQ has a barrier, then 

(2.5) \imu{x) = 

It is therefore important to understand which domains in the Euclidean plan satisfy the 
hypothesis of Theorem 12.31 Indeed, general sufficient conditions can be described in terms 
of local properties of the boundary (see for instance [58l Proposition 5.13]). 

Theorem 2.6 (Lebesgue). The Dirichlet boundary value problem H2.1\) is solvable for ar¬ 
bitrary assigned continuous boundary value values if every component of the complement of 
the domain consists of more than a single point. 


For the applications of this paper, the following corollary is essential. 
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Corollary 2.7. Let Q be a Jordan domain, then the Dirichlet boundary problem 112.1\) is 
solvable in for arbitrary continuous boundary values. 

The (strong) maximum principle (see for instance, [16]) implies that a strong solution is 
unique. For VL = A & polygonal planar annulus, we can therefore make the following 

Definition 2.8. We call u G C‘^{A) fl C^{A) the strong solution of the Dirichlet boundary 
value problem of the Laplace equation, if 


(2.9) 


Am = 0 in .4., 

M = 1 on El, and m = 0 on E 2 . 


We end this subsection by recalling the following dehnition which is valid for any harmonic 
function. 


Definition 2.10 (A smooth harmonic conjugate ([Ml Chapter 1.9])). Let {xq, yo) be a point 
in A, and let {x, y) in A be an arbitrary different point. Let 7 be a simple, counter-clockwise 
oriented, piecewise smooth curve joining (xq, yo) to (x, y) in A. Let (3 be any simple, closed, 
counter-clockwise oriented, piecewise smooth curve in A whose winding number is equal to 
1 , and let s denote the arc-length parameter of these curves. 

A harmonic conjugate of u is dehned by 

r Qii 

(2.11) u*{x,y) = u*{xo,yo) + J 


where u*{xo,yo) is some arbitrary hxed constant, and 


( 2 . 12 ) 


period(M*) 



Remark 2.13. It is well known that a smooth harmonic conjugate u* is dehned up to a 
constant, i.e., an assigned value at a point in the annulus; furthermore, the function values 
at any point differ by integral multiples of its period (see for instance [Ml Chapter 1.9]). 


2.2. Discrete harmonic and asymptotically harmonic functions, and their conju¬ 
gates. We now turn to dehning a combinatorial function analogous to u*. We will start 
with some notation and dehnitions from the subject of discrete harmonic analysis that will 
be used throughout the rest of this paper (see for instance [S] or [Ml Section 1.1]). Let 
r = iy, E, c) be a planar finite network, that is, a planar, simple, and hnite connected graph 
with vertex set V and edge set E, where each edge (x, y) E E is assigned a conductance 
c{x,y) = c{y,x) > 0. Let ViV) denote the set of non-negative functions on V. Given 
E C V, we denote by E‘^ its complement in V. Set V{E) = {u E 'P(C) : S{u) C E}, where 

S{u) = {x E V : u{x) 7 ^ 0}. The set 6E = {x E E^ : {x,y) E E for some y E E} is 

called the vertex boundary of E. Let F = E U 5E, and let E = {{x,y) E E : x E E}. Let 

T{E) = [E, E, c) be the network such that c is the restriction of c to E. We write x ~ ?/ if 

(x, y) E E. The following operators are discrete analogues of classical notions in continuous 
potential theory (see for instance [M] and ini)- 
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Definition 2.14. Let u G V{F). Then for x G F, the function 

(2.15) Am(x) = E c{x,y) {u{x) -u{y)) 

y^x 

is called the Laplacian of u at x. For x G S{F), let {yi,y 2 , ■ ■ ■ ,ym} G F be its neighbors 
enumerated clockwise. The normal derivative of m at a point x G 5F with respect to a set 
F is defined by 

j 

(2.16) —(F)(x) = ^ c[x,y){u{x)-u{y)). 

y^x, yGF 

Finally, u G V{F) is called harmonic in F C Id if Au{x) = 0, for all x G F. 


In order to define the class of discrete, asymptotically harmonic functions, we will use 
notation and definitions from Appendix lA. II in particular, every triangulation in this section 
is assumed to be quasi-uniform and consisting solely of nonobtuse triangles (see condition 
(VO) in Appendix IATD . 

Given such a triangulation T of dl, following [SO], Chapter 2], for each one of its Voronoi 
cells which is centered at x,, we let 

/ y/2 

(2.17) Aj = max/(Fi,) and A = max A*, 

\j&Ni '' ’J' J iGJ 

where /(•) denotes Euclidean length. 


Definition 2.18. Let a G M be a positive constant. Let T be a triangulation of A, with 
Voronoi cells A function g : ^ R is said to be asymptotically harmonic of 

order a with respect to conductance constants {cjj = c(xj,Xj)}, if for alH G J we have 

(2.19) lA(xi)| = I c{ij){g{xj) -^(xi))| = 0(A"). 

j&Ni 


The following lemma provides an estimate for the integral of the normal derivative of g, 
which can be thought of as the discrete flux of g through the boundary of a Voronoi cell. 


Lemma 2.21 (Asymptotic flux estimate). Let g : —)■ R 6 e harmonic or asymptotically 

harmonic of order a, with respect to conductance constants {ctj = Then, for any 

homotopically trivial (in A), closed path 7 C which contains an integer number of 
Voronoi cells VLi, i E J, we have 


( 2 . 22 ) 


' XG 7 


dg 

dn 


X = 



if g is harmonic, 

if g is asymptotically harmonic or order a. 


Proof. Let 12m = be the maximal collection of control volumes enclosed in 7 , and let 

Em be those edges of that lies in the interior of the bounded region enclosed by 7 . The 
first Green identity (see for instance [5|, Proposition 5]) implies that for u,v E V{Qm), we 
have 


(2.23) 








f du 

Au{x)v{x)+ —{nm)iy)v{y) 


GG7 


dn 
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We now let u = 1 in the above equality, and obtain 


(2.24) 


"" r Bn 

0 = ^Au{xi)+ —{nm)iy). 
»•_ 1 JvG'y 


i=i 

It therefore follows, by the definition of the combinatorial laplacian, that 

m 

(2.25) 

i=l j£Ni 

and the assertions of the lemma readily follow. 


r Pi 

0 = ^ ^ c(i,j)(M(j) -«(*)) + / ^(^m)(x). 


The following Corollary immediately follows. 


□ 


Corollary 2.26 (Asymptotic path independence). Let 71 and 72 be two simple paths in 
C A joining two vertices vi,V 2 G such that the path 7 ^^ o 7 ^ is trivial in 7 ri(. 4 ,), 

and contains an integer number of control volumes hi*. Then 
(2.27) 


' xG'yi 


Bn 


.X 


'XG72 


Bg 

Bn 


[X = 


0 , 


if g is harmonic, 


0(A“), if g is asymptotically harmonic or order a. 


Let 5 ^ be a discrete harmonic or asymptotically harmonic function. Before defining the 
notion of a combinatorial conjugate function for g, we will need to define a special class of 
paths in Thereafter, by integrating a generalized version of the normal derivative of g, 
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along a path from this class, the combinatorial conjugate function of g will be dehned at the 
vertices of the Voronoi cells of T. 

We let A denote the union of all Voronoi cells of T, and make 

Definition 2.28 (Flux fellow paths). Let Uq be a hxed vertex in a Voroni cell of T and 
let uj be any vertex in A. Let 7 a = [cuq, ■ ■ ■ ,ujk = oj] be a simple, piecewise linear curve 
in Ad) joining ojq to u. For each [a;j,a;j+i], z = 0,..., /c — 1, let n* be the vertex in T® 
on the unique edge intersecting [cj*, cjj+i], and which is to the right of [a;j,a;j+i]. Then 
7r = ['^^ 0 , • • •, Ufc-i] C Td) will be called the flux fellow path of 7 a. 

Remark 2.29. The discussion preceding condition (VO) in Appendin lA. II grants us that 77 - is 
indeed a path in Td); we orient each edge in 7 a according to increasing order of its vertices. 



Figure 2.30. The path 7 a (solid red) and 7 r (solid brown). 

With a;o,cu, 7 A and 7 r defined as above, we extend 7 a in Ad) to a simple, piecewise linear 
and homotopically trivial, closed curve in A, in such a way that the piecewise linear disk 
bounded by it, and 77 - are disjoint. We let F{uo, co, 7 a, 7 r) denote such a disk. Note that 7 a 
and F(a;o, cj, 7 a, 7 r) are not uniquely determined; however, once 7 a is chosen, 77 - is uniquely 
determined. (In Figure [230l F(a;o, ca, 7 a, 7 r) is the union of three Voronoi cells.) 

We now need to extend the notion of the discrete normal derivative (see Equation (I2.16p in 
Dehnition l2.14p in order to treat the special case arising in the definition above. Specihcally, 
vertices in which the formally defined normal derivative needs to be computed do not belong 
to the vertex boundary of the domain, rather, they are close to it (in a combinatorial sense) 
as is the case with vertices that belong to 77 - and 7 a, respectively. In the definition that 
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will follow, we will abuse notation and keep the notation in (I2.16p for the extension. The 
extension needed for the applications of this paper is determined by the following 

Convention. The computation of the integrand in 

(2.31) f |^(F(a;o,a;,7A,7r))(i/),i/ e 7r, 

Jye^T 

will include only vertices that are adjacent to y along an edge which intersects 7 a. 


It is time to make the third, and final key definition of this section. 


Definition 2.32 (A combinatorial conjugate). Let T be a triangulation of A, and let A 
denote the union of all Voronoi cells of T. Let cuq be a hxed vertex in A and let cu be any 
vertex in A, let 7 a be a simple, piecewise linear curve in A joining ujq to a;, let 7r be the flux 
fellow path of 7 a, and let F(a;o, w, 7 a, 7r) be a disk as described above. 

(i) Let 5 ^ be a discrete harmonic, or a discrete asymptotically harmonic function of order 
a. Then, a combinatorial conjugate of g is dehned by 

(2.33) f{uj) = T{ujq) + f |^(F(a;o,a;,7A,7r))(2/) for every cu G A, 


where g*{u:Q) is some arbitrary, fixed constant. 

We divide each 2-cell in A^^^ into triangles with vertices in A*'^^ and disjoint interiors. We 
then extend g* affinely over edges in A*^^) and triangles in A*^^b (The extended function will 
also be called the combinatorial conjugate of g.) 

(ii) Let a\ be any simple, counter-clockwise oriented, closed curve in A whose winding 
number is equal to 1. The period of g* is dehned by 


(2.34) 


period(( 7 *) 


! i&OLj- 


dn 


( 0 - 


Remark 2.35. A notion of combinatorial period for g* was introduced in [381 Section 2], and 
combinatorial provisions analogous to those in Remark l2.13l hold for g*. In addition, it is also 
clear that the dehnition is independent of the choice of A(a;o, cu, 7 a, 7r)) and of the choice of 
7a due to Corollary 12.261 

Remark 2.36. In the dehnition of g*, the main new diherences compared to |3H1 Dehnition 
2.6] are: (i) it incorporates the Voroni cells, and (ii) the integration of the normal derivative 
of g is computed along hux fellow paths. These are needed since in the applications of 
this paper, the combinatorial function analogous to u is asymptotically harmonic^ and not 
harmonic as was the case in [38] . 

Remark 2.37. The search for discrete approximation of conformal maps has a long and rich 
history. We refer to [l5] and m Section 2] as excellent recent accounts. We should also 
mention that a search for a combinatorial Hodge star operator has recently gained much 
attention and is closely related to the construction of a harmonic conjugate function. We 
refer the reader to |10| and to m for further details and examples for such combinatorial 
operators. 
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3. Uniformization of a planar annulus 

In this section, we prove the main theorem of this paper, Theorem 13.121 and afterwards, 
we will indicate how the hypotheses “polygonal boundary” in the main theorem, can be 
relaxed to “continuous boundary”. 

3.1. Uniformization of a planar polygonal annnlus. In this subsection, we will continue 
to assume that ^ is a fixed polygonal annulus. We will construct a sequence of mappings, 
obtained via a refined sequence of quasi-uniform triangulations and conductance constants 
along edges according to fl3.13p . from the interior of A onto the interior of a concentric Eu¬ 
clidean annulus. The dimensions of the image annulus are determined (see Equation fl3.16p i 
by a specific Dirichelt boundary value problem. Theorem 13.121 demonstrates that the se¬ 
quence converges uniformly on compact subsets to a conformal homeomorphism. 

We keep the notation of the previous sections and Appendices (A)-(C). In particular, let A 
be endowed with a family of quasi-uniform triangulations {7^„} such that pn 0, as n oo. 
For each Tp^, let the corresponding family of Voronoi cells be denoted by In 

addition to requiring that each triangulation is quasi-uniform, we will assume the existence 
of a constant tq such that for all pn < tq, the corresponding family of Voronoi cells 
satisfies the following assumptions: 

(VI): The number of essential neighbors of each Vi remains uniformly bounded, that 
is, 

max{cardAj} < m* for some m* G N; 

iaJ 

(V2): Each point Xjj = [xi^Xj] fl Tjj is the middle point of the segment Tjj. 

For W <Z A and two points a,b ^ W, we let D{a,b) denote the pseudo-distance on W 
defined as the infimum of the Euclidean lengths of curves in W that join a to b. We then 
define the intrinsic diameter of hi by 

(3.1) IDiam(iy) = sup{U(a, b) \ a, 6 G W}. 


In Lemma [3751 we will need to consider Dirichlet boundary value problems with prescribed 
Poisson data and non-homogeneous boundary condition as described in the following 


Definition 3.2. Let h : dA ^ M be a continuous function and let h be a C^{A)r\C^{AUdA) 
extension of h. A function u in is a (strong) solution of a Dirchlet boundary 

value problem for the Poisson equation and non-homogeneous boundary data, if for / = —Ah 
and for g = h\dA, we have 


(3.3) 


Am = / in and 
u = g, on dA. 


Remark 3.4. The existence and uniqueness of a strong solution to fl3.3l) follows from Corol¬ 
lary [2771 by setting u = u — h, where u is the strong solution of fl2.ip . 

Several ideas from the proof of the following lemma which provides an estimate for the 
discrete flux of a solution of fl3.3p along the full boundary of one Voronoi cell, will be essential 
in the proof of Theorem 13.121 We will apply a weaker version of the following 
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Lemma 3.5 ([30l Lemma 2.63]). Assume that conditions (VI) and (V2) hold. Let the so¬ 
lution of as. A) belong to C"^(r2). Then there exists some constant c = c(u,fl) such that 


(3.6) 


E 


- u{Xi)) 



< cAf, for all i E J, 


where A, is defined by 

Remark 3.7. Assumptions (VI) and (V2) are critical in the proof of this Lemma. In the 
statement of this Lemma, 12 is the domain on which fl3.3p is initially dehned. However, 
the proof remains valid even if the regularity assumption is only assumed to hold for any 
close, proper subset of 12 with thick enough neighborhood; that is, if the subset and its 
neighborhood are still contained in 12. This weaker assumption will be used in the proof of 
Theorem 13.121 below, where we will also show how to choose a thick neighborhood. 


Let us now restrict our attention to the particular type of boundary value problem fl2.9p 
that will be considered in Theorem 13.121 Let h be the continuous function dehned on dA by 
setting 

(3.8) h\Ei = 1 and h\E 2 = 0. 

Let h G C^{A) n C^{dA) be an extension of h to the interior of A. Let us denote the 
projection of h on 12„ by n„(h), that is, we hrst dehne 

(3.9) n„(h)(a;) = h(a;), for every x G 
and then extend affinely over triangles. 

Let n be a solution of fl2.9p . Recall that u = u + h, where the trace of h G L'^{A) is equal 
to the trace of u on dA, and u is the solution of the homogeneous Dirichlet boundary value 
problem for Poisson’s equation 


(3.10) 


Am = —Ah, in A 
u = 0, on dA. 


Definition 3.11. Let m and = Up^ be the solutions of flA.ip . presented as in flA.22p with 
/ = —Ah, respectively. 


With the notation above, we now turn to the main theorem of this paper. In the statement 
below, quasi-uniform triangulations are the subject of Dehnition IA.71 Dehnition lA.lOl and 
Equation flA.28p explain the terms appearing in the dehnition of the conductances. Equa¬ 
tion fl3.13l) . and exploits the connection to Stephenson’s conjecture (Conjecture 10.21) . In 
order to ease the notation, we will not distinguish between a map dehned on the 0-skeleton 
of a triangulation and the affine extension of the map on edges and triangles. Finally, for 
every n, let 7;^ = 
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Theorem 3.12. Let {%i} be a sequence of quasi-uniform triangulations of A of mesh size 
Pn ^ 0 as n ^ oo, and let the corresponding family of Voronoi cells of each Tn be denoted by 
Assume in addition that {Tn} satisfies conditions (1^1) and (V2). Let the conductance 
of each edge e ^ T, T ^ Tn^'^ be defined by 

(3.13) Cn{e) = 


ij,n 


d. 


i],n 


Let u and h be given in (12.91) and (13.8p . respectively, and define (see Definition LS.11\) ) 
(3.14) gn = Un + 'n.n{h). 

Then, as n ^ oo the following assertions hold: 

( 1 ) 11 '*^ ~ hnWiooiA) 0 . 

(2) On each proper, compact subset of A, the pn’s are asymptotically harmonic of order 
a = 3. 

(3) Let g* denote the combinatorial conjugate of pn, and let <pn be the sequence of discrete 
mappings defined by extending affinely the discrete sequence of mappings given by 

2-ji 

(pniu:) = exp {dniuj) + ^ e 

over fin- Then the sequence {(fn} converges uniformly on compact subsets of A to 
a conformal homeomorphism, denoted by onto the interior of the concentric 
Euclidean annulus £_a, whose inner and outer radii are given by 


(3.16) 


(3.16) 


^ {l,exp(- 


27r 


where u* and period{u* 


'period(M 
are given in Definition \2.1(X 


)}’ 


Remark 3.17. Following this procedure for every p > 0, we have turned Tp into a hnite 
electrical network, as predicted in Stephenson’s conjecture (see Conjecture 10.2p . where the 
induced potential function dehned on it (see [52]) satishes the system of equations described 
by flA.31D . We remark that since for each p > 0, the values of Up at the boundary vertices 
dT^ C dft are given, there is no need to specify conductance constants for edges that are 
contained in dfl] or one can choose arbitrary values. 


Proof. Setting / = —Ah in Theorem IA.231 we know that for the approximation of u by 
(see Equation fl3.10l) h the following estimate holds. There exist constants C = C{A) and a 
real number s = s(.4.) e (0,1) such that 

(3.18) \\u-Un\\L^{A) < C'Pnlog(l/Pn)||Ah||i^(^). 

In particular, as p„ —?■ 0 we have that 

(3.19) 11“ ~ “n||l,oo(A) 0. 

Since both h and n„(h) are continuous in A, we also have that 

(3.20) ||h - n„(h)||L^(^)-)■ 0. 
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Hence, we are now able to show that the gnS comprise our desired approximations to m - 
the strong solution of the smooth Dirichlet problem for the Laplace equation fl2.9p . Precisely, 
by recalling Remark 13.41 we observe that 

\\U- (Unih) +Un)\\L^{A) = \\u - h + fl - {Un + Iln{h))\\ L^(A) 


(3.21) 


U + h- {Un + 'n.n{h))\\L^(A) 


||(m Un) (h n„(h)) 


< \\u - Unh^iA) + \\h - n„(h)||L^(^). 

Therefore, the hrst assertion of the Theorem is proved by applying Equations (I3.19p and 

(1^ . 

We will continue the proof by showing that for n large enough, each gn is asymptotically 
harmonic of the same order. To this end, recall that, by dehnition, for every vertex Vi E 
we have 

(3.22) Iln{h){vi) = h{vi). 

It is clear that h is a solution of fl3.3p with g = h\dA and with / = Ah. Therefore, for 
each n, applying Equation fl3.6p with hn = n„(h) and / = Ah, together with taking p = 
and Up = Un in Equation flA.31D . imply that each g^ is discrete asymptotically harmonic of 
order a = 3. By construction, it also holds that gn\'Tn^'’ fl Ei = 1 and gn\'T^ fl ^2 = 0, for all 

n. 

We now turn to prove the uniform convergence of the ^*’s, over compact subsets of A, to 
u*- the harmonic conjugate of u. Indeed, let C ^ be a compact annulus with smooth 
boundary which is concentric with A, where e = dist(9^, 9^*^) is small. Let us choose n 
large enough (i.e., small enough) so that there exists a triangulation Tp^ E {Tn} satisfying 
the following. 

(V3): A subset of its associated volume elements is contained in A^, and 

the combinatorial one vertex neighborhood of this subset, when considered in 
is also contained in A"^. 

Following Dehnition 12.281 and the discussion proceeding it, we choose an index i G J', as 
well as one of the vertices of which will be denoted by ojq. Let oj be any vertex in 

Let = [ojq,uji, ... ,uJk-i,oj = ojk] be a (piecewise linear) simple, counter¬ 
clockwise oriented path in the one skeleton of Ap^ = joining cjq to cj. We will also 

dehne u*(uo) = 0 . 

It then follows from Dehnition 12.101 and Remark 12.131 that the smooth harmonic conjugate 
function u* satishes 

f du 

(3.23) u*{oj) = / ^ds. 

We now follow the notation in Dehnition 12.281 and we let 77 ^^ denote the hux fellow 
path of 7 A^^[a;o,a;]. Let us write 7 rp„ = K",..., and set TnM = 0. Hence, by 
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Definition I2.32[ the asymptotic combinatorial conjugate of gn = gp„ is defined at u by 
(3.24) gJM = [ 


We now turn to estimate 
(3.25) |t<*(ij)-9;(ij)| 


'7Ap„ 


du , 
on 


'®G7r„, 


dgn 

dn 


(X 


In fact, we are interested in this difference as n —?> oo. Hence, by the extension of the 
definition of the combinatorial normal derivative (see Equation fl2.16p and Equation fl2.3ip ). 
and the first assertion of the theorem, it is sufficient to replace gn, in any term appearing in 
the second integrand in fl3.25l) . with the restriction of u to the vertices in 

It is classical that u G and therefore Equation (5.11) in [30] shows that for 

Xi = Xi^p„,Xj = Xj^p^, Tij = and with A* = (see fl2.17p ). there exists a positive 

constant cq = co(m,AI^) so that 


(3.26) 


J-(u(xj) -u(xi)) 

Ui,j 


du 

dni,j 



< CqA^, i G J',j G iVj- 


Equation (5.15) in [30] shows that for each i G J', the following continuous linear functional 
Ti (which can actually be defined for any function in C^(r2i,p„)) 


(3.27) 




j&Ni 


the following estimate holds 

(3.28) 

where ci = ci{u,Qi^pJ. 

It therefore follows that 

(3.29) 

f du 


\Tiu\ < ciA- 


'7A„„ 


dn 


ds — 




dgn 

dn 


[X 


k „ 

■r, J \u> 


< 


< 


E 

i=0 


IDiam(Al'^ 

A 


du , 

du 

{Xij)) 

dni,i+i 

ll 

du , 

du 


dni,i+i 

UTI i,j 


■C2 (m, Al^)A^ 


< C‘i{u,A^)\ ^A- 

< csA^. 


It is evident that as n ^ 0 we have that A —?■ 0. 

We now study several choices that are to be made as we let n ^ cx3, and as a result the 
triangulation is changing. Let m > n, and let Tm = Tp^ be the corresponding triangulation 
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which we may again assume satisfies condition (V3) with the appropriate indices changes. 
Following the paragraph preceding (V3), we let go and q be the analogous choices for ojq and 
uj, respectively. Let a be any smooth path in joining uq to go, we then define 

f du 

(3.30) r = r(cr) = j -^ ds , 
and accordingly set 

(3.31) g*^{qo) = r. 

Then, it readily follows that an estimate analogous to Equation fl3.29p . holds for the 
difference 

(3.32) \n'(q)-S’M\- 

Hence, we conclude that 


(3.33) 


lim 

n^oo 


sup 


xGUiO.. 


(0) 


u*{x) - g*^{x)\ 


0 , 


and as n —)■ cx), with i G J'{pn) comprises a dense subset of A^. Therefore, applying 

the uniform continuity of the g*’s and u* in A^ and letting e —?■ 0, we obtain uniform 
convergence of the g*’s to u* over compact subsets of A. 


We also need to prove that 


(3.34) period(g*) —)■ period(M*). 

To this end, let us choose a point Pq in and let {3 and period(M*) be given according to 
Definition 12.101 Furthermore, let ujn G be chosen so that Un ^ Pq- Let 'jj-n a closed 
curve in An\ based at oon according to which period(g*) is computed. 

Since u* is continuous in A'^, we have that 


(3.35) u*{ujn) —)■ u*{Po) as n —)■ cx). 

By applying now Equation fl3.32p with we can conclude that fl3.34p holds. 
It now follows that the (p^s converge uniformly on compact subsets of A to 

277 - 

(3.36) $^(z) = exp (m(z) + zu*{z ))). 


We end the proof by recalling a classical result (see for instance [20l Section 7] or [581 
Theorem 4.3]) which asserts that <F is a conformal homeomorphism between the interiors of 
A and respectively. 

□ 

The case of continuous boundary. In this paragraph, we will briefly indicate why the 
boundary regularity assumption in Theorem 13.121 can be relaxed. Assume that .4, is a planar 
annulus, where dA is a union of disjoint, non-degenerate, continuous Jordan curves. 
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Definition 3.37 ([HI 1.6.7]). We will say that a sequence of planar annuli IZj G TZ, j = 
1,2,..., with {rj, r|} as the components of their complements converges from the inside to 
an annulus 7Z with {T^i, 7 ^. 2 } as components of its complement, if for every e > 0 there exists 
such that for n> every point of {T^j)i=i ^2 lies within a spherical distance less than e of 
the set {TZi,TZ 2 )i=i, 2 - 

A classical construction due to Kellogg im Chapter XL 14] grants us an existence of a 
nested sequence of annuli, {Ai}, where for all i > 0, {Ai} C A, the boundary of Ai is 
polygonal, and the sequence converges to A from the inside. Furthermore, since each At is 
made of a lattice of squares, it is easy to construct a sequence of qausi-uniform triangula¬ 
tion each Ai which satisfy assumptions (V0)-(V3). Thus, A is presented as an increasing 
union of open subsets, the interiors of the Ai, and each conformal embedding can be 
approximated according to Theorem 13.121 It follows that (up to normalization of the maps 
<h^J, a subsequence of the will converge uniformly on compact subsets of A, to its 

uniformizing map (see for instance m Lemma 2.2] or [U Page 223 ]). Hence, we have the 
following 

Corollary 3.38. With the additional approximation processes described in the paragraph 
above, we may assume in Theorem \3.1S\ that dA is continuous. 

4. The simply connected case 

In this section we will affirm Stephenson’s conjecture 1 Conjecture 10.2p which was originally 
stated for the case of a bounded, simply connected, planar domain. Our point of departure 
is Theorem 13.121 whose notation will be closely followed. The proof of this case entails 
successively applying this theorem to an increasing sequence of annuli, a known modifica¬ 
tion of Koebe’s compactness theorem, Riemann’s removable singularity theorem, a lemma 
concerning the monotonicity of periods, and a basic covering property of planar Riemann 
surfaces. 

In the following, we will let 
(4.1) a{z) = -, 

be the standard inversion of C; it is well known that a is conformal. We can now turn to 

Theorem 4.2. Let Q be a simply connected domain, embedded in C and bounded by a closed, 
polygonal curve T; let po E Ll be a fixed point. Let {Hn} C Ll be a nested seguence of disjoint, 
polygonal, Jordan disks with polygonal boundaries {0n} such that the disks converge to po, 
that is, 

(4.3) 111 D H 2 ... D flfc ..., 

(4.4) meshifiln) ^0 as n ^ 00 , and 

(4.5) Po = n„Lln- 

For each n, let An = .4.„(L2, 0„) be the polygonal annulus defined by Q\Qn with dAn = 
T U 0„, endowed with a seguence of guasi-uniform triangulations {Tm,An}m=i’ 
all m = m{An) large enough, Tm,An satisfies the hypotheses of Theorem VJ.ltl Let 

(4.6) : An —)■ £-n 










20 


SA’AR HERSONSKY 


be the sequence of conformal homeomorphisms constructed according to Equation / (,?. onto 
the interior of concentric Euclidean annuli £n, whose inner and outer radii are given by, 
respectively 

2 % 

(4.7) ^^)}. 

where m* is the (smooth) harmonic conjugate ofun, the solution of the boundary value prob¬ 
lem h2.9i) defined on An- 

Then, a normalized subsequence of {a o converges uniformly on compact subsets of 
Q\Pq to a holomorphic homeomorphism ^ from Et\pQ onto D \ 0. Eurthermore, ^ can be 
extended to be holomorphic overQ. 

Proof. By construction, the {An} is a strictly increasing sequence, that is, 

(4.8) C ^2 C ... C A ... 

which all share F = dQ as their outer boundary component, and with hi \ {po} being their 
union. The following lemma is needed in order to understand a monotonicity property of 
the sequence {An}- 

Lemma 4.9. The sequence {period{u}()} is strictly decreasing. 

Proof. By Green’s theorem, for all n > 1 we have that, 

(4.10) f I'Vunl'^dx[ AunUndx = [ ^^ds. 

JAn JAn JdAn 

However, by the dehnition of period(M„), and since Un is the solution of the boundary value 
problem fl2.9p dehned on An, for all n > 1, we have that 

(4.11) / \Vun\^dx = period(M„). 

It is clear that for all n > 1, can be extended to be zero on An+i \ An+i to a piecewise 
smooth function on An+i having the same boundary values as those of Un+i- The assertion of 
the lemma now follows by the well-known characterization of Un+i as the unique minimizer 
of the Dirichlet integral over An+i- 

Lemma IT^ 

It follows from Equation (14.Ih . Equation (14.7p and the Lemma, that the sequence {An = 
a(£n)} consists of planar, concentric. Euclidean annuli, such that the inner and outer radii 
of each An are given by 

27r 

(4.12) = {l/exp(—^).l}. 

respectively; where the sequence {r 2 ,n} is strictly increasing. Note that all the H„’s share 
§1 = (9D as their outer boundary component, 

(4.13) AiCA2C...CAk..., 
and the sequence {An} exhausts D \ 0. 
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Figure 4.14. The evolution of 11. 


Pick zo e Ai and a local complex parameter at zq. For all n > 1, we now apply a 
normalization by post composing a o vl>„ with a conformal embedding /„ : —)■ C so that 
the composed maps 

(4.15) = /„ o a o ; X -t C 
satisfy 

(4.16) T„(zo) = ^0 and T'^{zo) = 1. 

Note that the image of each An is still a concentric Euclidean annulus, yet the sequence 
{T„(^„)} is not (generically) concentric. Nevertheless, it follows from a modihcation of 
Koebe’s compactness theorem (see for instance [T^ Proposition 7.5]) and a Cantor diago- 
nalization process, that a subsequence of the {T„} converges uniformly on compact subsets 
of 11 \po) to a conformal, univalent mapping 

(4.17) T:ll\po^C, 

which is obviously not constant. It is also evident that T is bounded, and therefore, by Rie- 
mann’s removable singularity mapping theorem, can be extended to a conformal, univalent, 
embedding from H. Hence, the extended map must be equal to the Riemann mapping with 
the same normalization. This ends the proof of the Theorem. 

□ 
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Remark 4.18. Following the rationale preceding Corollary 13.381 we can assume that F is 
continuous. 


We conclude the main body of this paper with a few comments. 

1. Disk packing and quasi-uniform triangulations. It is well known (see for instance 


Section 5]) that a sequence of disk packings satisfying some minor conditions, induce 
(as explained in IA.8P a sequence of quasi-uniform triangulations, that will in addition satisfy 
assumptions (VO) and (VI). However, assumption (V2) will not always be satished; it will 
be satisfied (for instance) for sub-packings of scaled copies of the infinite hexagonal disk 
packing (which were the subject of Thurston’s original conjecture). Recall that assumption 
(V2) was used in the proof of Theorem 13.121 only in the part addressing the convergence of 
the We will relax this assumption in 


2. The case of higher connectivity. As mentioned in the introduction, Stephenson’s 
original conjecture can be formulated for any finitely connected, Jordan domain. However, 
several issues need to be addressed before an appropriate statement can be made. For 
instance, the existence of singular points and level curves for smooth harmonic functions 
solving a Dirichlet problem (analogous to the one in Theorem I3.12p on such domains needs 
to be addressed. We will treat this issue and others in 


Appendix A. preparatory facts 

A.l. The finite volume method. In this section, we will assume that H is a fixed, 
bounded, m-connected m > 2), polygonal domain in The vertex centered finite volume 
method (FVM) is a powerful discretization scheme, aimed at constructing and presenting 
approximations of solutions of partial differential equations in the form of algebraic set of 
equations, where the unknowns are placed at the vertices of a discrete grid. The phrase “fi¬ 
nite volume” refers to control volumes associated with a chosen neighborhood of each vertex. 
This method turns out to be useful in the case of elliptic problems of diffusion type, such as 


(A.l) V ■ (KVu) = / in H, / G L^(H), 

with prescribed boundary conditions on dQ. Here, the diffusion coefficient K{x) is a 2 x 2 
symmetric matrix with components in L°°{Q), which is uniformly positive definite in H. 

By integrating over H and employing Gauss divergence theorem, it follows that the equation 
above can be written as 

(A.2) - [ (KVu) ■nds= [ f dx, 

Jdn Jn 

where h denotes the unit outwards normal vector of dkl, ds denotes arc length along dVt, 
and dx is the (standard) area measure in the plane. 

A common theme in finite volume methods is to approximate the solution of Equa¬ 
tion flA.2p by replacing it with integrating Equation flA.ip over control volumes, while taking 
into account the boundary conditions. Let us now describe an important special case of a 
boundary value problem which will be essential to the applications of this paper. We start 
with 
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Definition A.3. A triangulation T of is a set of triangles Tj, i = 1,... ,?7, such that the 
following hold. 

n 

(A.4) Tj, and Tj n Tj = 0, a vertex or one common edge, for all i ^ j. 

i=l 

The following quantity is associated with a fixed triangulation. 

Definition A.5. Let T be a triangulation on hi, the mesh size of T is equal to 

(A.6) sup (i(T), 

Ter 

where d{T) denotes the diameter of T (i.e., the length of its largest edge). Henceforth, Tp 
will denote a triangulation of H of mesh size equal to p. 

In order to apply the machinery of numerical approximation of elliptic boundary value 
problems, one needs to avoid situations where triangles in (any) Tp become fiat as p ^ 0. 
To this end, we let cr(T) denote the diameter of the largest circle that can be inscribed in a 
triangle T. We now define the geometric property of the class of triangulations that will be 
used throughout this paper. 

Definition A.7. A family of triangulations {Tp} of H is called quasi-uniform (or quasi¬ 
regular) when p —)■ 0, if there exists a positive constant r such that 

(A.8) [ I < r, for all T G 7),, and for all p. 

a{T) 

The induced control volumes, or the Voronoi cells used in this paper, which are associated 
with T are defined as follows. For each triangle T E T, let ct denote the circumcenter of T, 
which by definition is the intersection point of the perpendicular bisectors of the edges. We 
join ct' to Ct by a segment [ct', ct] whenever T and T' share an edge. This procedure divides 
each (interior) triangle T into three quadrilaterals and induces a new decomposition of H. 
The control volume of a vertex n G T is defined to be the star of v in this decomposition. 

Henceforth in this paper, we will assume that 

(VO): every quasi-uniform triangulation T under consideration is quasi-uniform and 
consists exclusively of nonobtuse triangles. 

It then follows (see for instance [3l Theorem 6.5]) that T is a Delaunay triangulation, i.e., 
no point in the vertex set of T lies inside the circumcircle of any triangle in T, and the 
corresponding Voronoi diagrams can be constructed by means of the perpendicular bisectors 
of the triangles’ edges. 

Let {Tp}p>Q be a family of triangulations of H, and let the set of vertices of Tp be enu¬ 
merated by {ui,^ 2 ) • • • Wm(p)}- denote the set of vertices of T E Tp, and let 

Vp{Tp) denote the set of interior vertices of Vp{Tp) = UrGrpb]o(T). We now turn into inte¬ 
grating Equation (IA.2p over each control volume Vt^.,i = 1,... ,M{p). For each vertex Vi, 
i = 1,..., M{p), let Aj denote the set of indices of neighboring vertices of n*. Each is an 
open, simply connected, and polygonally bounded set. Its boundary, dfly^, consists of finitely 
many (straight) line segments Fj^j = fl j = 1,... ,ni, where Ui is the number of 
vertices adjacent to Vi, along which the normal n\Tij = fiij is constant. 





24 


SA’AR HERSONSKY 




Figure A.9. A triangle and a control volume. 


Definition A.10. Let mjj denote the length of Fjj, and let dij = \vi — Vj\, where |, | denotes 
the Euclidean distance in 


Note that by the discussion preceding Definition lA.lll the length of Fij is equal to |ct—ct' |, 
where T and T' are the (only) two triangles that Tij intersects. The union is called 

the Voronoi diagram of Tp. 

We now consider the homogeneous Dirichlet boundary condition defined in 12 by letting 
K be the identity matrix in Equation (lA.lj) . and 0 = 0 as the boundary condition. Given a 
triangulation T of 12, the finite dimensional space defined below is used. 


Definition A.11. 


(A.12) Vo,r = {n : G ^ M |n e C'(G), v\T e Vi{T) for all T e T and n = 0 on 
where Vi (T) denotes the space of linear polynomials in two variables over T. 


One important feature of Vo,r is that it is a linear subspace of a certain Sobolev space in 
which the weak solution of the above boundary value problem is constructed. Let us recall 
the following definition. 

Definition A.13. The Sobolev space iL^’^(12) is the subset of L^(12) defined by 
(A.14) = {ve L\Q) I d^v, dyV e L\Q)}, 

where dxV, dyV denote the distributional derivatives of v in the x and the y directions, re¬ 
spectively. The integration is with respect to the standard Lebesgue measure in the plane 
which will be denoted by dx. 

For G iL^’^(12), one defines the scalar product and an associated norm, respectively, 
by 

(A.15) {u,v)i^ 2 = / {uv + Vu-Vv), |m|i 2 = (m,m)i ,2 

Jo. 

where Vn = {d^v, dyv), and the scalar product is the Euclidean one in It is well known 
that iL^’^(12) equipped with this scalar product is a Hilbert space. Finally, let be 

defined as the closure of (^“(G) in with respect to this norm. Equipped with this 

scalar product iLQ’^(r2) is a Hilbert space as well. It is a fact that Vo,r is a linear subspace 
of Hq (12), whose dimension is equal to the cardinality of the 0-skeleton of T. 

The first step in finite volume methods amounts to finding Up G "^o,Tp such that 
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(A.16) / Vwp • n ds = / / dx, for all v G Vp{Tp), 

J dQv J 

where h denotes the outer unit normal to dVly. In order to ease the notation, we will suppress 
the index p when a hxed triangulation Tp is under consideration. Next, one studies the 
convergence in (various norms) of the tip’s to ti, as the mesh of each triangulation converges 
to zero. This analysis is crucial to the applications of this paper and in the next section we 
will recall an approximation result from this theory. 

A.2. Piecewise linear approximations of the solution. In this section, we will follow 

closely the exposition in [HI Sections 1-3] and apply a special type of finite element method 
in the case of homogenous Dirichlet boundary value problem, with the goal of writing it in 
a variational form. 

Let 

(A.17) Fp = {?7 G L 2 (fl)) I r]\VLy is constant for v G VpiTp), and r]\VLy = 0 if n G cAl}. 

The following is the Petrov-Galerkin formulation of the hnite volume method of Equation 
(IA.2|) . Find Up G Vo, 7 ^ such that 


(A.18) ap{up,p)-=- r]{v) f Wup ■ h ds = f fp dx, for all r; G Fp. 

vGV0{Tp) 

Let Xv be the characteristic function of Qy, and let Ip : G(f2) —)■ Fp be the interpolation 
operator dehned by 

(A.19) Ipfi= Y (t>{v)Xv 

Recalling Dehnition lA.lll we note that Equation flA.16p can be transformed to its varia¬ 
tional form 
(A.20) 

j{up,Ip(j)) = - Y / Vup-hds= Y / / dx, for all (/) G Vo.Tp- 


OCr 


veVO{Tp) vevO(Tp) 

Note that for every f ^ Lp and G Vq^t^, the following holds 


-T/ \ do do„ 


(A.21) 

•el'.ra V<EV«(T,) 

Let US denote the right-hand side of Equation flA.2ip by (/, Ipfi). Hence, Equation flA.20p 
can be written in the form 


(A.22) ap(Mp, Ipfi) = (/, Ipfi) for all G Vo,rp- 

The existence of the solution Up follows from the coerciveness of the inner-product ap for 
p, small enough (see for instance [13] L 
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For the applications of this paper, it is necessary to consider non-convex polygonal do¬ 
mains. We will construct sequences of such domains with uniform upper and lower bounds 
on their largest and smallest angles, respectively, in order to approximate (in a sense that 
we will explain later) a given Jordan domain. In particular, due to the presence of corner 
singularities of vertex angles that are bigger than vr, the recent L^o error analysis of the finite 
volume element which is needed for our applications is quite subtle (see for instance Ha)- 

We need to introduce some notation before stating the main result which we will be 
using. Let hi be a bounded, (possibly) non-convex, polygonal domain. It is well known, 
that if / G Lp(r2), 1 < p < oo, then the solution u of the boundary value problem in 
flA.lj) is not always in H^{Q) (see for instance [29l Section 2]). However, it turns out that 
u always belongs to a fractional order Sobolev space for some 0 < s < 1; where s 

is effectively determined by the maximal interior angle of H and p (see [TH Section 2] for 
precise definitions). 

The following foundational result that was obtained by Chatzipantelidis and Lazarov is 
the main numerical approximation result which will be used in this paper. 

Theorem A.23 ([TH Theorem 4.8]). Let u and Up he the solutions of HA . 1\} and M. J61) . 
respectively, with f G Lp,p > 1. Then there exists a constant C, independent of p, such that 


(A.24) 


M-Mplkoe < C'P*log(l/p)ll/IUp- 


A.3. Stephenson’s conductance constants - a finite volume method perspective. 

In this section, we will express each Up as the potential of an electrical network. To this 
end, two steps are necessary. First, we will present the form Op in terms that relate to the 
geometry of the given triangulation Tp. Second, we will turn the triangulation into a finite 
network. The second step amounts to assigning conductance constants along the edges of Tp. 
The values of the conductance constants will be extracted from the first step. 

Thus, we will see that Up can be presented as the solution of a system of finitely many 
linear equations. Finally, from this presentation, we will define a quantity, the discrete flux 
of Up, which will be important for the applications of this paper (see [3HIEZ] and [Ml Section 
1] for a background on hnite networks and discrete flux in the above setting). 

In [1] , it was shown that under the assumption that K is the identity matrix the following 


holds 



(A.25) 


Hence, Equation flA.22D takes the form 

(A.26) f Vup -Vf dx= (/, Ipf) for all f G Vo,rp- 

Jn 

(It is well known that for every p > 0, Vq^t^ is a hnite dimensional vector space which is 
spanned by {4>i}p (the nodal basis). Therefore, it is sufhcient to check the equation holds 
for any two elements in the nodal basis.) 

We keep the notation as in the discussion preceding Figure IA.91 The following Lemma 
and its corollary allow us to turn each Tp into a hnite network. 
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Lemma A.27 ([H Lemma 6.8]). Let Tp be a triangulation of Vt, and consider its corre¬ 
sponding Voronoi diagram. Then, for an arbitrary triangle T ETp with vertices Vi,Vj{i ^ j), 
the following relation holds 


(A.28) 


• V0J dx = 


m 




di 


JT ^13 

where mjj is the length of the segment ofTij that intersects T. 


A computation then shows that 


Corollary A.29 ([H Corollary 6.9]). 


Under the assumptions of Lemma\A.2'?\ we have 


(A.30) 



Hence, by letting the index i range over the indices of the interior vertices (i.e, those that 
are in VpifTp)), flA.26p turns into the following system of linear eguations 


(A.31) 




/ dx. 


Remark A.32. Note that when / = 0, Up is a discrete harmonic function on with the 

conductance constant ^ for the edge joining Vi to n,-. 

dij j 
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